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ABSTRACT 

We report the detection of high energy 7-ray emission from the young and energetic pul- 
sar PSR B1509— 58 and its pulsar wind nebula (PWN) in the composite supernova remnant 
SNR G320.4— 1.2 (aka MSH 15—52). Using 1 year of survey data with the Fermi-Large Area Telescope 
(LAT), we detected pulsations from PSR B1509— 58 up to 1 GeV and extended 7-ray emission above 
1 GeV spatially coincident with the PWN. The pulsar light curve presents two peaks offset from the 
radio peak by phases 0.96 ± 0.01 and 0.33 ± 0.02. New constraining upper limits on the pulsar emis- 
sion are derived below 1 GeV and confirm a severe spectral break at a few tens of MeV. The nebular 
spectrum in the 1 - 100 GeV energy range is well described by a power-law with a spectral index of 
(1.57 ± 0.17 ± 0.13) and a flux above 1 GeV of (2.91 ± 0.79 ± 1.35) xlO'^ cm'^ s"!. The first 
errors represent the statistical errors on the fit parameters, while the second ones are the systematic 
uncertainties. The LAT spectrum of the nebula connects nicely with Cherenkov observations, and 
indicates a spectral break between GeV and TeV energies. 

Subject headings: gamma rays: observations - ISM: individual (G320.4— 1.2, MSH 15— 52) - pulsars: 
individual (PSR B1509-58, PSR J1513-5908) 
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1. INTRODUCTION 

Pulsars and pulsar wind nebulae (PWNe) are be- 
lieved to be sources o f cosmic-ray electrons (e.g., 
IKennel fc COTOTdtil flOSl IGaensler fc Slmg 120061) . Al- 
though hadronic 7-ray emission from TeV emitting 
PWNe has been suggeste d by many authors (e.g., 
iBednarek fc Bartos ik 2003; Horns! [20061) . most evidence 
indicates that 7-rays are generated via inverse Comp- 
ton scattering of electrons accelerated in pulsar magne- 
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tospheres and at pulsar wind termination shocks. 

The composite supe rnova remnant (SNR ) G320.4— 1.2 
(aka MSH 15-52; iCasweU et all Il98l is usually 
associated with the rotation-powered radio pulsar 
PSR B1509-58 (aka PSR J1513-5908). The 150 ms 
rotation period was disc overed by the Einstein satellite 
(jSeward fc Hardenlll982D and soon thereafter confirmed 
in the radio domain |Manchester et al]|1982f ). With a 
large period derivative (1.5 x 10~^^ss~^), this pulsar is 
one of the youngest and most energetic pulsars known in 
the Galaxy with a characteristic age of 1700 years and a 
spin-down power E of 1.8 X 10^7 erg s'^ The inferred 
surface magnetic field is 1.5 x 10^^ G derived under the as- 
sumption of a dipolar magnetic field. The measurement 
of the pulsar braking index n = 2.839 shows that thi s 
assumption is reasonable (e.g.. [Livingstone et al.|[2005l) . 
Therefore, the high magnetic field is not much below the 
quantum-critical magnetic field of 4.413 x 10^'^ G, the 
domain of the so-called high-B-field pulsars and mag- 
nctars. The distance is estimate d at 5.2 ifc 1.4 kpc us - 
ing HI absorption measurements ([Gaensler et al.lFl999l) . 
This is consistent with the value of 4.2 ifc 0.6 kpc derive d 
from the dispersion measure (iCordes &: Laziol [200I . 
PSR B1509— 58 has been studied by all major X-ray 
and 7-ray observatories yielding a broad-band spectral 
energy distribution and pulse profiles as a function of 
energy. Its detect ion by COMPTEL (0.75 - 30 MeV, 
iKuiper et al.l |1999() and non-detection with the Ener- 
getic Gamma-Ray Experiment Telescope (EGRET) in 
the 30 MeV - 30 GeV energy range, both aboard the 
Compton Gamma Ray Observatory (CGRO), indicate 
an abrupt spectral break between 10 and 30 MeV. This 
break is well below the break energies of most 7-ray pul- 
sars detected by Fermi which are typically around a few 
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GeV (|Abdo et al.ll2010aD . More recently, the detection 
of pulsed 7-rays from PSR B1509 — 58 at 4 a level abov e 
100 MeV was reported by AGILE (|Pellizzoni et al.ll2009[ ). 

Einstein X-ray observations of MSH 15—52 revealed 
an elonga ted non-t hermal source centered on the pul- 
sar (S eward fc Harden.1982 ). later confirmed by ROSAT 
and interpreted as a pulsar wind n ebula powered by 
PSR B1509-58 (|Trussom et al.lfT996l) . This PWN, com- 
posed of arcs and b ipolar jets, is especially bright and 
extended in X- rays (Tamura et al." 1996; iGaensler et alj 
120021: iForot et al.. .2006 : Yatsu ct al., 20091). and at very - 
high energies (ISako et a l. 2000; lAharonian et al.l 120051 : 



iNakamori et all l2008f ) . The dimensions o f the PWN as 
observed bv ROSAT (Tlhissoni et al.l ll996h and H.E.S.S. 
(|Aharonian et al.ll2005D are 10' x 6' and 6.4' x 2.3' re- 
spectively. The multi- wavelength emission of the PWN 
in MSH 15—52 can be accounted for by synchrotron radi- 
ation from electrons within the PWN and inverse Comp- 
ton scattering on soft photons such as the cosmic mi- 
crowave background (CMB), the infrared (IR) aiid the 
optical interstellar rad iation field (|Aharonian et al.]|2005l : 
INakamori et al.|[2008[ ). This leptonic model requires a 
broken power-law spectrum for the electrons. However, 
large uncertainties on the break energy remain due to 
the lack of observations at corresponding wavelengths, 
namely IR or optical for the synchrotron radiation and 
GeV 7-rays for the inverse Compton component. The 
observations performed now with Fermi help constrain 
the electron spectrum, in particular the break energy, of 
the PWN of MSH 15—52 and provide new elements to 
the discussion on the energetics. 

Successfully launched in June 2008, the Large Area 
Telescope (LAT) onboard the Fermi Gamma-ray Space 
Telescope covers the 20 MeV - 300 GeV energy range. 
With its improved performance compared to its prede- 
cessor EGRET, it offers the opportunity to search for 
high energy pulsations of PSR B1509— 58, and, measur- 
ing the spectrum of the PWN in MSH 15—52, to better 
constrain the emission models in pulsar winds. The re- 
sults of 1-year observations of the pulsar PSR B1509— 58 
and its PWN are reported in the following sections. 

2. RADIO TIMING OBSERVATIONS 

With its large spin-down power, the pulsar 
PSR B1509— 58 is a good candidate for 7-ray de- 
tection an d is monitored b y the LAT pulsar timing 
campaign fS mith et aLll2008D coordinated among Fermi, 
radio and X-ray telescopes. 

The ephemeris of the pulsar PSR B1509— 58 used in 
the analysis of the Fermi-LAT data is obtained using 
observations at 1.4 G Hz madc^ with t he 64-m P ar kes ra - 
dio telescope ()ManchesterH2008l : IWeltevrede et alll2009[ ). 
A total of 42 time of arrivals (TOAs) were recorded be- 
tween 2007 April 30 and 2009 August 29. Radio observa- 
tions simultaneous with the 7-ray data make it possible 
to correct for the large drift in phase caused by timing 
noise when constructing the 7-ray light curves. 

The TEMP02 timing package (|Hobbs et al.ll2006f ) is 
then used to build the timing solution from the 42 
TOAs. We fit the TOAs to the pulsar rotation frequency 
and its first two derivatives. The fit further includes 3 
harmonically-related sinusoids, using the "FITWAVES" 
option in the TEMP02 package, to flatten the timing 
noise. We used the value of DM = (252.5 ± 0.3) cm^^ pc 



for the dispersion measure (|Hobbs et al.ll200l . given m 
the ATNF Pulsar Catalogue The post-fit rms is 875 
/xs, or 0.6% of the pulsar phase. This timing solution 
is used in the temporal analysis described in detail in 
Section liXn 

3. LAT DESCRIPTION AND DATA SELECTION 

The LAT is a high energy photon telescope sensitive 
to 7-rays with energies from below 20 MeV to more 
than 300 GeV, that detects photons through pair conver- 
sion. The photon incident direction is derived by track- 
ing the electron-positron pair in a high-resolution con- 
verter tracker and the energy is measured with a CsI(Tl) 
crystal calorimeter. An anticoincidence d etector identi- 
fies t he background of charged particles (jAtwood et al.l 
l2009f ). In comparison to its predecessor EGRET, the 
LAT has a larger effective area (~ 8000 cm^ on-axis), a 
broader field of view (■^ 2.4 sr) and a superior angular 
resolution (^ 0.6° 68% containment at 1 GeV for events 
converting in the front section of the tracker). 

The analyses reported here are performed on 374 days 
of data taken in survey mode (2008 August 04 - 2009 
August 13). Events from the "Diffuse" class are selected, 
i.e. the highest quality photo n data, having the m ost 
stringent background rejection (|Atwood et al.ll2009| ). In 
addition, we exclude events with zenith angles greater 
than 105° to avoid contamination by 7-rays produced by 
cosmic-ray interactions in the Earth's atmosphere and 
periods corresponding to a rocking angle (i.e. the angle 
between the viewing direction of the LAT and the zenith) 
larger than 43°. 

4. RESULTS 

4.1. Gamma-ray emission from PSR B1509—58 
4.1.1. Light curves 
We selected photons with an angle 9 < max(5.12° x 

energy of the pho- 



'lOOMeV 

ton, from the radio pulsar position, R.A. ^ (228.48175 ± 
0.00038)°, D ec. = (-59.13583 ± 0.00028)° in J2000 
(|Kaspi et al.l Il994f ) The energy-dependence of the in- 
tegration radius is a satisfactory approximation of the 
shape of the LAT Point Spread Function (PSF), espe- 
cially at low energies. 

Photons in this energy- dependent region are then 
phase-folded using the timing solution described in Sec- 
tion [2] The resulting 7-ray light curve for energies higher 
than 30 MeV, is presented in Figure[TJ We have a total of 
28966 7-rays in the circular region of energy-dependent 
radius, among which are 1267 ± 515 pulsed photons after 
background subtraction. The radio profile (red dashed 
line) obtained from the 42 radio observations included 
in the timing solution used for our analysis is overlaid 
in Figure [T] for comparison. In this analysis, phase is 
defined as the maximum of the main radio peak observed 
at 1.4 GHz. 

In Figure (TJ two peaks PI and P2 can be ob- 
served at phases 0.96 ± 0.01 and 0.33 ± 0.02 respec- 
tively. The uncertainty in the extrapolation of the ra- 
di o pulse arrival time to infinite frequency, as defined 
in iManchester and Tavloil (|19770 . is 0.004 in phase and 

http;/ /www. atnf.csiro.au/rescarch/pulsar/psrcat/ 
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Fig . 1. — Light curve of the pulsar PSR B1509— 58 above 30 MeV within an energy-dependent circular region, as described in Section 
14.1.11 The light curve profile is binned to 1/30 of pulsar phase. The radio profile (red dashed line) is overlaid in arbitrary units. The main 
peak of the radio pulse seen at 1.4 GHz is at phase 0. Two cycles are shown. 



can be neglected. Hence, the peaks are separated by 
= 0.37 ± 0.02. PI and P2 are symmetric and can be 
well modeled by Lorentzian functions of half-widths of 
0.22 ± 0.11 and 0.05 ± 0.03 respectively. We also notice 
that the first 7-ray peak leads the radio main pulse by 
phase 0.04 ± 0.01, as shown in Figure [I] which corre- 
sponds to a delay of (6 ± 2) ms. 

Table [1] pres e nts th e value of the H-test as defined in 
Ide Jager et"all (|1989( ) and obtained in the 30 MeV - 100 
GeV, 30 MeV - 100 MeV, 100 MeV - 300 MeV, 300 
MeV - 1 GeV and 1 GeV - 100 GeV energy bands us- 
ing the energy-dependent circular radius defined above. 
The corresponding light curves are presented in Figure [5] 
along with the light curve measure d by GOMPTEL i n 
the 0.75 - 30 MeV energy range (|Kuiper et all 1 1 9991 ) . 
Within the error bars, the peak positions remain stable 
with energy. From Figure [5] and Table [1] we notice that 
no significant pulsation can be detected above 1 GeV. 
Using 374 days of data in survey mode, a H-test value of 
31.34 is obtained in the 30 MeV - 100 GeV energy range, 
corresponding to a significance of 4.51 cr. 

4.1.2. Spectral analysis of PSR B1509-58 

A spectral analysis of the pulsar is performed in the 
100 MeV - 1 Ge V energy range usi ng a maximum- 
likelihood method (jMattox et al.lll9^6l ) implemented in 
the Fermi SSG science tools as the gtlike code. This 
tool fits a source model to the data along with models 
for the instrumental, extragalactic and Galactic back- 



Energy band 


H-test 


Significance 


(GeV) 




(^) 


0.03 - 100 


31.34 


4.51 


0.03 - 0.1 


15.42 


3.07 


0.1 - 0.3 


15.60 


3.09 


0.3 - 1.0 


4.66 


1.42 


1.0 - 100 


0.06 


0.03 



TABLE 1 

Results of the periodicity test applied to PSR B150 9— 58 

USING THE energy-dependent REGION DEFINED IN SECTION |4. 1. 1! 



grounds. In the following spectral analysis, the Galactic 
diffuse emission is modeled using the ring-hybrid model 
glLiem_v02.fit. The instrumental background and the ex- 
tragalactic radiation are described by a single isotropic 
component with a spectral shape described by the tab- 
ulated model isotropic-iem_v02.txt. These models and 
their detailed description are released by the LAT Col- 
laboration Sources near the pulsar PSR B1509— 58 
found above the background with a statistical signifi- 
cance la rger than 5 a are ex tracted from the source list 
given m lAbdo et al.l (|2010bD . and are taken into account 
in this study. We use P6_V3 post-launch instrument re- 
sponse functions that take into account pile-up and ac- 
cidental coincidence effects in the detector subsystems 
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Fig. 2. — Light curves of the pulsar PSR B15 09— 58 in different energy bands within a circular region of energy-dependent radius. From 
bottom to the top: COMPTEL (0.75 - 30 MeV MKuiper et al.llT999l ') and LAT profiles in 30 MeV-300 GeV, 30 MeV-100 MeV, 100 MeV-300 
MeV, 300 MeV— 1 GeV, 1 GeV-300 GeV energy bands are presented. Two cycles are shown. 



65 

Despite the detection of pulsations, no significant 7- 
ray emission can be observed at the position of the 
pulsar using 374 days of LAT data. Indeed, most of 
the weak signal detected on PSR B1509— 58 is observed 
at low energy (below 300 MeV) where the LAT an- 
gular resolution is large in comparison to the distance 
that separates our source of interest from the bright 
pulsar PSR J1509-5850 (less than 0.8°). This ren- 



ders the spectral analysis extremely complex. There- 
fore, only 2 (T upper limits can be derived and are pre- 
sented in Figure |31 As an attempt to evaluate the flux 
of PSR B1509— 58, pulsed excess counts were derived 
from the light curves (presented in Figure [5]) in the 
100 - 300 MeV and 300 MeV - 1 GeV energy bands, 
as well as the corresponding effective area. The rapid in- 
crease of the effective area in the 100 - 300 MeV energy 
range (Atwood ct al. 200^ makes the analysis highly de- 
pendent of the assumed spectral shape and yields very 
65 See |http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/cia¥ggij^fi-i[§gEpiaa^L^"g-Ig^^ Only the 

for more details ^ 
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result obtained in th e 300 MeV ~ 1 GeV energy range 
(jAtwood et al.l[2009l ). where the effective area is more 
stable, is represented in Figure [31 The upper limit de- 
rived with gtlike in this energy band is consistent with the 
integrated flux estimated from the number of pulsed pho- 
tons. The overall spectrum in the 1 keV - 1 GeV energy 
range indicates a very low cut-off o r break energ y in th e 
pulsar spectrum, as suggested by iKuiper et al.l ()1999t ). 
The implications of such upper limits on the emission 
models in pulsar magnetospheres are further discussed 
in Section [5Tl 



4.2. High energy analysis of the PWN in MSH 15—52 
4.2.1. Morphology 

Figure m presents the smoothed counts maps of the re- 
gion around MSH 15—52 in Galactic coordinates above 
1 GeV (top panel) and 10 GeV (bottom panel) and 
binned in square pixels o f side-length 0.05°. The H.E.S.S. 
(jAharonian et al.l l2005l ) contours have been overlaid in 
black for comparison. At low energies, the emission 
is essentially do minated by the bri ght nearby pulsar 
PSR J1509-5850 (jAbdo et al.ll20T0a[) . marked by a blue 
star, whereas the significant 7-ray emission above 10 GeV 
is spatially coincident with the nebula in MSH 15—52. 
The position of the pulsar PSR B1509— 58 is also marked 
with a blue star. 

An analysis tool. Sourcelike, developed by the LAT 
Collaboration allows us to estimate the position and the 
size of the source, assuming a spatial and spectral model 
for the diffuse emission and different morphologies: a 
point source, a Gaussian shape and a uniform disk. In 
this method, the likelihood is iterated to the data set as- 
suming spatial source models, taking into account nearby 
sour ces, Galact ic diffuse and isotropic components in the 
fits (lAbdo et a l. 2GlGc), as described in Section l4X2l 

Morphological studies are performed above 6.4 GeV, 
assuming the three spatial hypotheses mentioned above. 
The choice of this energy is motivated by the better an- 
gular resolution, the non-contamination from the Galac- 
tic diffuse emission, and nearby bright sources such as 
the pulsar PSR J1509— 5850. The positions, extensions 
and the corresponding errors, as well as the Test Statis- 
tics (TS) obtained for each hypothesis are summarized in 
Tableia The TS values are obtained as TS = 2(Li -Lq), 
where and Li are the values of the log-likelihood ob- 
tained by null hypothesis and each source hypothesis, 
respectively. In view of the errors of localization, the fit 
positions are compatible with each other, and bests fits 
are obtained either with a uniform disk (TS of 69.4) or 
a Gaussian distribution (TS of 67.6). 

The differences in Test Statistics between the two ex- 
tended shapes and the point source hypotheses: TSe^t 
— 24.7 and TSext — 22.9 for the uniform disk of exten- 
sion (T ^ 0.25° and the Gaussian shape of extension a ^ 
0.15° respectively with relation to a point source, indi- 
cate a significant extension for the LAT source spatially 
coincident with the PWN in MSH 15-52. The conver- 
sion of the extensions obtained for the uniform disk and 
Gaussian distribution into a rms value gives results con- 
sistent with each other. In the following analyses, the 
spectral results obtained for the extended scenarios are 
presented. For comparison, the results derived assuming 
a point source hypothesis are also quoted. 



4.2.2. Spectral analysis 

The following spectral analyses are performed using gt- 
like. Sourcelike, described in Section [4.2.11 is also run as 
a cross-check and gives compatible results. The mod- 
els used to describe the Galactic, extragalactic and in- 
strumental components are mentioned in Section 14.1.21 
Sources near the PWN with a s tatistical significanc e 
larger than 5 a are extracted from lAbdo et al.l (|2010bl ) , 
and are taken into account in this study. 

Since the pulsar PSR B1509— 58 is seen only below 
1 GeV, the 7-ray photons in the 1 - 300 GeV energy 
range, in a 20° x 20° square centered on the the pulsar 
radio position and coming from the entire pulse phase 
interval are selected. 

The spectral analysis is performed above 1 GeV assum- 
ing different morphologies for the source: a point source, 
a Gaussian shape and a uniform disk using the positions 
and extensions summarized in Table |2l A spectral fit 
using the H.E.S.S. image of the PWN (|Aharonian et al.l 
[2005li as a template is also performed. In this energy 
range, for all spatial distributions, the LAT spectrum 
can be well modeled by a simple power-law: 
dN / \-r 



cm^^s^^MeV"^ 



dE " VlGeVy " ™' 

where F is the photon index of the spectrum. The inte- 
grated fluxes and spectral indices of the source obtained 
for different spatial hypotheses are summarized in Ta- 
ble [3l The best fit is obtained for a uniform disk hy- 
pothesis, favored with respect to the Gaussian, H.E.S.S. 
template and point source morphologies with differences 
in TS of 3.9, 13.7 and 32.9 respectively. The LAT spec- 
trum for a disk hypothesis is well described by a power- 
law with a spectral index of (1.57 ± 0.17 ± 0.13) and a 
flux above 1 GeV of (2.91 ± 0.79 ± 1.35) x lO"'^ cm^ s'^. 
The flrst errors represent the statistical error on the fit 
parameters, while the second ones are the systematic un- 
certainties. 

Three different uncertainties can affect the LAT spec - 
trum estimation, as described in lAbdo et"al] (|2010d|) . 
The first one is due to the uncertainty in the Galac- 
tic diffuse emission since MSH 15—52 is located close 
to the Galactic plane. Different versions of the Galactic 
diffuse emission generated by GALPROP (|Strong et al.l 
were used to estimate this error. The differ- 
ence with the best diffuse model is found to be less than 
6%. This implies systematic errors on the fluxes of 33% 
above 1 GeV. The second systematic is related to the 
morphology of the LAT source. As described in Ta- 
ble m the Gaussian, disk hypotheses and the H.E.S.S. 
template match the gamma-ray morphology quite well. 
The fact that we can not decide which one is better 
adapted induces an additional systematical error on the 
flux of the order of 24% above 1 GeV. The third system- 
atic is produced by the uncertainties in the LAT instru- 
ment response functions (IRFs). We bracket the energy- 
dependent effective area with envelopes above and below 
the nominal curves by linearly connecting differences of 
(10%, 5%, 20%) at log(£;/l MeV) of (2, 2.75, 4) respec- 
tively, which yields additional errors on the flux and spec- 
tral index. 

The spectral energy distribution (SED) of the PWN 
in the case of a uniform disk hypothesis is presented in 
Figure [5l The Fermi-LAT spectral points were obtained 
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Fig. 3. — Spectrum of PSR B1509— 58 from soft X-rays to 7-rays. 2 a upper limits derived with gtlike from LAT observations are 
represented by red stars. The green diamond represents the integrated flux derived from the light curves and estimates of the effective 
area and exposure in the 300 MeV - 1 GeV energy range, as described in section 14.1.21 The filled squares are the COMPTEL flux points 
derived from the excess counts in the 0.15 - 0.65 phase range while the open square represents the 10 - 30 Me'V flux in the 0.15 - 0.65 phase 
interval above the spatially determined background. The blue circles are the 2 a upper limits for the total fluxes obtained by EGRET. 
The different lines represent the best spectral fit measured by ASCA (short-dashed line, 0.7-10 keV), Ginga (dot-long-dashed line, 2-60 
keV), CGfiO-OSSE (long-dashed line, 50-750 keV), WELCOME (dotted line, 94-240 keV) and RXTE (dot-short-dashed line, 2-250 keV) 
IIKuiper et al.|[T999l , and references therein). 



Spatial model 


Galactic 


Galactic 


Error 


Radius 


Test Statistics 




longitude 


latitude 


(degrees) 


(degrees) 






(degrees) 


(degrees) 








Point Source 


320.288 


-1.209 


0.028 




44.7 


Gaussian 


320.275 


-1.266 


0.051 


0.146 ± 0.023 


67.6 


Uniform disk 


320.270 


-1.271 


0.061 


0.249 ± 0.047 


69.4 



TABLE 2 

Position of the centroid, extension and significance of the PWN in MSH 15—52 obtained with Sourcelike applied to the LAT 

DATA. 



by dividing the 1 - 100 GeV range into 6 logarithmically- 
spaced energy bins and performing a maximum like- 
lihood spectral analysis in each interval, assuming a 
power-law shape for the source. These points provide a 
mo del- independent maximum likelihood spectrum, and 
are overlaid with the fitted model over the total energy 
range (black line). 2 a upper limits are derived in energy 
bands where the significance level of the signal is lower 
than 3 a. 



Spatial model 


Flux above 1 Ge'V 
(lO-^cm-^s-i) 


Spectral index 


Point Source 


2.00 ± 0.76 


1.57 ± 0.24 


Gaussian 


3.01 ± 0.81 


1.58 ± 0.17 


Uniform disk 


2.91 ± 0.79 


1.57 ± 0.17 


H.E.S.S. 


2.22 ± 0.77 


1.52 ± 0.21 


TABLE 3 



Spectral parameters of the PWN obtained with gtlike for 
different spatial models. Statistical errors only are 

QUOTED. 



5. DISCUSSION 



5.1. High energy emission from PSR B1509—58 

The improved performance of the Fermi-JjAT com- 
pared to its predecessor EGRET allows the first detection 
of pulsations from PSR B1509-58 up to 1 GeV with a 
light curve presenting two peaks at phases 0.96 ± 0.01 
and 0.33 ± 0.02 as seen in Figure [T] The alignment 
of the broad peak (P2) at phase <j> — 0.33 is in gen- 
eral agreement with measurements of the phase of the 
main pulse by other high-energy exp eriments, e.g. 0.38 ± 
0.03 with COMPTEL (10 - 30 MeV.'Kui per et al.|[T999t) . 
0.32 ± 0.02 with CGRO-BATSF. and CGiZO-OSSE (20 - 
500 keV, lUlmer et al.l Il993fl and 0.35 with AGILE 
(0.1 - 30 GeV. lPellizzoni et al.ll2009D . In X-rays, the peak 
shifts towards phase 0.24 - 0.27 (e.g. iKawai et al.l 119911: 
iRots et al.iri998l: ICusumano et all 1200 It) . We do not de- 
tect the suggested narro w pulsed coinponeii t located at 
phase 0.85 reported bv lKuiper et al.l ()1999D using com- 
bined COMPTEL and EGRET data in the 10 - 100 MeV 
energy range. However, we observe a pulse component at 
phase 0.96 ± 0.01. Considering the marginal significance 
of the COMPTEL/EGRET pulse, these peaks might be 
one single component. AGILE also observe d a possible 
peak at phase ^ 0.85 (|Pellizzoni et al.ll2009[ ). which ap- 
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Fig. 4. — Smoothed counts maps in arbitrary units of the region around MSH 15—52 above 1 GeV (top) and 10 GeV (bottom) 
in Galactic coordinates, binned in square pixels of side- length 0.05°. The maps are smoothed with a Gaussian of cr = 0.15°. The 
H.E.S.S. (Aharonian et al. 2005) contours of the PWN in MSH 15—52 are overlaid in black for comparison. The positions of the pulsars 
PSR J1509-5850 and PSR B1509-58 are marked by blue stars. 
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pears to be broader than the COMPTEL pulse. The de- 
tection of the narrow pulse only above 10 MeV strongly 
suggests a harder spectrum for this component than for 
the broad one. 

The Fermi pulse shapes can be modeled with symmet- 
ric Lorentzian functions with half widths of 0.08 ± 0.06 
and 0.21 ± 0.09 for PI and P2, respectively. Below MeV 
energies, the broad component is asymmetric. Its pro- 
file can be described with 2 Gaussians com ponents (e.g. 
iKuiper et aLlll999t ICusumano et all 1200 It ). These com- 
ponents peak at phases ~ 0.25 and ~ 0.39 with widths 
of 0.056 and 0.129, respectively. Above a few tens of 
MeV, it seems that the first of the two components is 
no longer contributing to the pulse profile, causiii g the 
apparent 'shift', as observed bv lKuiper et al.l (|1999t ): the 
COMPTEL broad peak can be described with one single 
component. Within errors, the widths of the broad peak 
measured by COMPTEL and Fermi are in agreement. 

The pulse component observed at phase 0.96 ± 0.01 
leads the 1.4 GHz radio pulse by 0.04 ± 0.01 in phase. 
This feature is quite remarkable as so far all radio- 
loud pulsars excluding millisecond pulsars detected by 
Fermi present phas e lags with respect to the radio pulse 
(jAbdo et al.ll2010al) . However, the first 7-ray peak of the 
Crab pulsar precedes the main radio peak, b ut the high- 
energy peak lags the radio precursor (e.g. lAbdo etakl 
I2010d ). In the case of PSR B1509-58 a we ak com- 
ponent was reported by ICrawford et al.l ()2001[ ) at 1351 
MHz, which precedes the main pulse by ~ 0.14 in phase. 
This might be a precursor as is seen in the Crab, but this 
has to be confirmed before we draw conclusions based on 
this possible feature. 

Before Fermi launch, one of the major open questions 
dealt with the zone of gamma-ray production in pul- 
sar magnetospheres. There are two classes of models 
that differ by the location of the emission region. First, 
polar-cap (PC) models place th e emission near the mag- 
netic poles of the neutron star (|Daughertv and Harding! 
[1996). The second class is fornied by the outer-ga p 
(OG) models (e.g. ICheng et al.l [19861: IRomanil [1993), 
where the emission extends between the null charge sur- 
face and the light cylin der, and by the two-po le caus- 
tic (TPC) models (e.g. iDvks and RudakI |2003| ) which 
migh t be realized in slot gap (SG ) acceleration models 
(e.g. iMuslimov and Harding! [2004[ ) . where the emission 
takes place between the neutron star surface and the light 
cylinder along the last open field line. After one year of 
Fermi observations, the high-altitude models seem to be 
favored, even if these models do not work for all pulsars. 
From previous high-energy observations, it was not ob- 
vious beforehand whether the high-energy emission from 
PSR B1509— 58 could be described by surface or outer- 
gap emission. 

The 7-ray pulse profil es can be used to co nstrain the 
geometry of the pulsar. IWatters et al.l (|2009l ) simulated 
a population of young spin-down-powered pulsars for 
vacuum-dipole magnetospheres. The peak separation of 
the 7-ray pulse profile presented in Figure [T] and the ra- 
dio lag can yield constraints on the viewing angle C and 
the magnetic inclination a. 

Using the heuristic law for the 7-ray luminosity 



(jWatters et al.ll2009D : 



« T]E w C X 



E 



1/2 



2033 ergs 



xlO^^ergs"^ (2) 



with C a slowly varying function of the order of unity, the 
7-ray efficiency t] is estimated at 0.007. Assuming a 7- 
ray peak separation of 0.37 ± 0.02 and a 7-ray efficiency 
of 77 = 0.01, the clos est value to the real 7 -ray efficiency 
that can be found in IWatters et al.l (|2009f) . yields a tight 
constraint for C in the range 64° — 70° and a values of 
45° — 65° in the framework of the OG model. The pulse 
profile can also be explained with the TPC model assum- 
ing C of 50° — 65° and a of 45° — 65°. However, the radio 
delay with respect to the first 7-ray peak could not be 
explained neither by the TPC nor the OG models. 

Constraints on a and C are also possi- 
ble with radio-polar iz ation measurements (e.g. 
Radhakrishnan fc Cookd 119691 : iLvne and Manchester 



.19881) . Polarization measurements by ICrawford et al . 
(}2001l) show highly linearly polarized signals of 97% 
and 94% at 660 MHz and 1351 MHz, respectively. 
However, the position angle shows only a shallow swing, 
which suggests a large magnetic-pole impact angle 
[3 = ( — a. Magnetic inclination angles larger than 60° 
are excluded, in agreement with the maximum values 
of /? derived using the OG and TPC models of 25° and 
20°, respectively. 

From the spectral analysis, the stringent upper limits 
of the pulsed spectrum measured by Fermi-h AT and pre- 
sented in Figure [HI confir m the spectral break between 10 
and 30 MeV suggested bv lKuiper et al.l ([1999! ). However, 
the new 2 a upper limits obtained with Fermi give more 
stringent constraints than EGRET's one. The emission 
m odels can now b e tested against these constraints. 

iRomani! (|1996l) modeled the high-energy emission 
based on curvature radiation-reaction-limited charges in 
the outer magnetosphere and argued that, for high- 
field pulsars such as PSR B1509— 58, synchrotron flux 
will dominate the emission in the 100 keV - 10 MeV 
band, and more specifically the GeV band curvature 
component. The reason for this is that high altitude 
two-photon pair creation in collisions between X-rays 
and hard gamma-rays is prolific, enhanced by aberra- 
tion effects. This process will curtail super-GeV flux 
while permitting synchrotron emission at lower ener- 
gies by electrons with ~ 0.3 - 10 GeV energies. In 
this way, the radiative power is transferred in cascading 
from the super- GeV band to a < lOMeV synchrotron 
window, i.e. generating a spectrum peaking below the 
Fermi-LAT energy range. Based on the three dimen- 
sio nal outer-magnetosp h ere model of pu ls ars p roposed 
by ICheng et al.! ([20001 ) . iZhang fc Chei^ ([200I calcu- 
lated the light curves and spectra of PSR B1509— 58 as- 
suming a = 65° and ( = 75°. Their light curve presents 
a single broad peak comparable with the RXTE pulse 
profile, though a bit too narrow ([Rots et al.!li998[ ). The 
resulting spectrum is characterised by a simple power- 
law from the soft X-ray band to a few hundred keV, with 
a photon index of 1.5 and a cut-off below 1 MeV. This 
overall shape agrees well with the multi- wavelength data 
available even though the proposed cut-off lies at an en- 
ergy s lightly smaller than that proposed bv lKuiper et al.l 
([19991) which is consistent with the new LAT data. Fur- 
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thermore, their light curve presents a single broad peak 
while a second peak is observed with the LAT data above 
3 MeV. 

iHarding et al.l ()1997t ) argued that the polar cap model 
could explain the low spectral cut-off observed by CGRO 
for PSR B1509— 58. Curvature emission at low al- 
titudes would naturally appear at energies above 10 
GeV, which was obviously not seen by EGRET and 
cannot be discerned in the LAT data presented here. 
Spectral attenuation by magnetic pair creation 7_B — > 
e+e~ would be extremely effective at en e rgies above 
100 MeV in the scenario of IHarding eiTall (fl99l . Us- 
ing Eq. [1] of iBaringI (|2004[ ). the maximum photon en- 
ergy consistent with magnetic pair production trans- 
parency at altitude r along the last open field line is 
-Emax - 1.76(Bi2)-iFi/2 (r/i?,)^/2 GcV, for a surface 
polar field strength of Bq ~ 10^^i3i2G and stellar radius 
i?, - lO^cm. Bo - 1.5 X 10^3 G and P = 0.15 s then 
set E'niax ^ 45(r/i?,)'^/2 MeV for the energy of the -yB 
pair creation turnover. Even at the stellar surface, this 
estimate is too high to accommodate the EGRET upper 
limits and COMPTEL data downturn at around 10 MeV 
as seen in Figure |3l 

However, as emphasized bv lHarding et al.l()1997D . mag- 
netic photon splitting, a quantum-electrodynamic pro- 
cess important only for magnetic fields approaching the 
quantum critical value Bcr = 4.413 x 10^^ G, can atten- 
uate 7-rays e niitted near the surfac e of strongly magne- 
tized pulsars. IHarding et al.l ()1997D showed that photon 
splitting will be important for 7-ray pulsars having a sur- 
face magnetic field larger than O.^B^r, where the splitting 
attenuation lengths and escape energies become compa- 
rable to or less than those for pair production. Specif- 
ically, they demonstrated that attenuation due photon 
splitting would reduce i?max to nicely accommodate the 
CGRO observations, but only if the emission was pre- 
dominantly at r i?, and its co-latitude was consis- 
tent with a standard polar cap in PSR B1509-58. The 
new Fermi-hAT data with the combined soft X-ray to 
soft 7-ray (COMPTEL) spectral points confirm this pic- 
ture: the polar-cap model is spectroscopically viable for 
this pulsar, but subject to the strong constraint of emis- 
sion at the magnetic co-la t itude of the rim, i.e. ~ 2° as 
proposed bv iKuiper et al.l (|1999D . Higher altitudes and 
accompanying larger co-latitudes will push -Emax to ener- 
gies above 50 MeV. Accordingly, the LAT suggestion of 
modest pulsations up to energies almost as high as 1 GeV 
indicates that some portion of this emission might em- 
anate from altitudes well above the stellar surface, where 
photon splitting will play a minimal role. 

Although there are only upper limits on the pulsed 
spectrum in the LAT energy range at this point, the emis- 
sion from both components of the light curve seems to ex- 
tend to 1 GeV. The broad peak at phase 0.33 is consistent 
with outer magnetosphere geometry and is also roughly 
in phase with the COMPTEL peak. If it is assumed that 
the radio peak arises at small magnetic colatitudes, then 
the radio/soft 7-ray phase separation suggests that both 
the Fermi and COMPTEL emission components at this 
phase originate in the outer magnetosphere, but given 
the sharp cut-off just above COMPTEL energies they 
must have different mechanisms. The narrow Fermi peak 
at phase 0.96 just leading the radio peak has no counter- 



part at COMPTEL energies but given its extension to 1 
GeV and the magnetic pair and photon splitting attenua- 
tion limits discussed above, this peak must also originate 
in the outer magnetosphere. However, its phase location 
is not easily explained by current outer magnetosphere 
gap models. High altitude pair-starved polar cap emis- 
sion has been showi i to produce a single peak just lead- 
ing the radio pulse (jVenter. Harding fc Guillemot) l2009l) 
and could potentially explain this component of the light 
curve at a similar a and C range. However, explaining 
both Fermi light curve peaks would require both pair- 
starved and non-pair-starved (gap) models to co-exist. 

One might expect the high magnetic field of 
PSR B1509— 58 to be the main reason for its unique be- 
havior. This seems not to be the case if one compares 
PSR B1509— 58 with the other pulsar Fermi has detected 
with an inferred surface magnetic field above 10^'^ G : the 
nearby (1.4 kpc) PSR J0007-h7303 (|Abdo et al.l [200I 
I2010al) has a period of 314 ms and a period deriva- 
tive of 3.61 X 10^'^ ss~^. Its inferred surface mag- 
netic is 1.1 X 10^^ G. The similar magnetic fields of 
PSR J0007-K7303 and of PSR B1509-58 seem not be 
the criterion for similar behavior as they appear to be 
very different. Indeed, PSR J0007-1-7303 belongs to the 
top ten brightest 7-ray pulsars and has a hard power-law 
spectrum with photon index F = 1.38 with a spectral 
cutoff energy at 4.6 GeV. Its spin-down power (4.5 x 10^^ 
ergs~^) is almost two orders of magnitude less than 
PSR B1509-58. Also PSR J0007-t-7303 shows no strong 
pulse d radio and X -ray emission like PSR B1509— 58 
tHalp ern et al.|[2"004( ) , is about ten times older and has a 
higher 7-ray efficiency than PSR B1509— 58. Finally, the 
pulse profiles and spectral behavior of PSR J0007-I-5303 
can be explained nicely with an outer-magnetosphere 
model unlike PSR B1509-58. 

5.2. Constraints on the emission models in the nebula 

High energy photons coming from pulsar magne- 
tospheres are usually expected to have a power-law 
spectrum with an exponential cut-off at a few GeV 
(|Abdo et al.ll2010a^ . iKuiper et all (|1999| ) suggested that 
PSR B1509— 58 has a spectral cut-off or break below 
100 MeV. Therefore, the absence of a pulsed signal above 
1 GeV, the spatial coincidence and the similar extension 
of the LAT source with the pulsar wind nebula as seen 
in X-rays and very high energy 7-rays, strongly suggest 
that the unpulsed 7-ray emission detected by the LAT 
above 1 GeV is dominated by the PWN. 

There are two possible interpretations for the origins 
of 7-ray photons from pulsar wind nebulae, i.e hadronic 
(from proton-proton interactions) or leptonic (via the in- 
verse Compton process). The multi- wavelength picture 
of MSH 15—52 is presented in Figure [6] using all avail- 
able d ata on MSH 15—52 as reported in Nakamori et al] 
()2008f ). The Fermi-L AT spectral points, obtained from 
the analysis described in Section 14.2.21 provide new con- 
straints on the m odel parameters . A s imple one-zone 
model described in lNakamori et aP (|2008[ ) can be used to 
reproduce the multi-wavelength spectrum of the PWN. 
We use the publicly distribu ted interstellar radiation 
field (ISRF) as described in iPorter et all (|2005i ) as tar- 

GALPROP model of cosmic- ray transp ort: 
|http: / / galprop.stanford.edu/web-galprop / galpr op-home.htmll 
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get photons for inverse Compton scattering. The ISRF 
spectra are modeled and given as a function of cylindri- 
cal coordinates in the Galaxy, for IR photons from inter- 
stellar dust grains, optical light from normal stars, and 
CMB. We do not consider the production of 7-rays via 
bremsstrahlung because of the low density of this regio n 
(up to ~ 0.4 cm"'^) as reported in lDubner et al.l ()2002t ). 
For simplicity, escape, energy and adiabatic losses as well 
as the time evolution of the magnetic field strength in the 
PWN are neglected, since the characteristi c age of the 
pulsar is quite young. As reported by .Nakamori et al.l 
HqOS), a single power-law electron spectrum does not 
reproduce the SED; hence the accumulated electron spec- 
trum used here follows a broken power-law with an ex- 
ponential cut-off: 



dE °^ 1 + {E/E^,)P2-P^ 



exp 



E 



E„ 



(3) 



where E^^^, i?br, Pi and p2 are the maximal energy, 
break energy and the indices of the electron spectrum 
respectively. The best fit yields the parameters listed in 
Table m and is overlaid in Figure IH] 



parameter 


value 


Pi 


1.5 


P2 


2.9 


Ebr(eV) 


4.6 X 10" 


-Bmax (eV).. 


2.6 X IQi-* 


B (mG) 


17 


(ergs) 


3.0 X 10** 


Wp (ergs) 


1.2 X 105i(1.0/n) 


Lay (erg/s)... 


3.9 X lO^'* 


Lie (org/s)... 


3.1 X 10^^ 



TABLE 4 

Parameters derived from the multi-wavelength spectral 
modeling 



The fitted mean magnetic field strength of 17 is 
identical with the valu e suggested by T eV observations 
(|Aharonian et al.l l2005t iNakamori et all [2008) and con- 
sistent with the lower limit of ^ 8 /iG obtained from 
X-ray data (jGaensler et al.l [20021 ) . We confirm that the 
7-ray emission from the PWN is dominated by the in- 
verse Compton scattering of the IR photons from in- 
terstellar dust grains with a radiation density fixed at 
1.4 eV cm""^ , which is the nominal value of the GAL- 
PROP ISRF. Here, the contribution of the optical pho- 
ton field is negligible because of the Klein-Nishina effect 
(jKlein fc N ishina 1929). 

The nature of the break energy Sbr remains unclear. A 
radiation cooling break which is dominated by the syn- 
chrotron loss in this leptonic scenario is a primary candi- 
date. With the assumed magnetic field {B = 17 /xG) and 
system age (r = 1.7 kyr), the bre ak energy can be c al- 
culated using standard formulae in iPacholczvkl (|197Cl[ ) as 
Q-Knilc' / aTT ~ 24 TeV leading to a break in the pho- 
ton spectrum at 36(S/1 iiG)(tIW yrs)-^ eV ~ 0.2 keV. 
Even considering the limit of the one-zone approximation 
and potential uncertainty in the parameters obtained, 
the fitted value of 460 GeV, which predicts a break in 
the photon spectrum at ~ 8 x 10^"^ cV, is much lower 
than the expected value. However, the photon index 



change AP — 0.7 is compatible with the prediction by 
the usual synchrotron cooling AP = 0.5. One should 
note that the synchrotron break in this pulsar wind neb- 
ula is expected to occur just below the X-ray band and 
no hints of a break have yet been observed. More likely, 
this break energy could be an intrinsic characteristic of 
the electron s pectru m injected in the PWN as suggested 
by Ide J aged (|2008D . though p2 = 2 is expected. An- 
other possibility would be that the spectral break of 
electrons may correspond to the energy scale where the 
electron acceleration mechan ism switches, for instance , 
from magnetic reconnection (jZenitani fc Hoshinol [200I 
to th e usual first-order F e rmi a cceleration. However, 
while IZenitani fc Hoshinol (|2001h indeed predict an in- 
dex of unity, the production of an index harder than 2 is 
explained with difficulty by Fermi acceleration. 

The total energy PSR B1509— 58 can supply to its 
PWN is strongly dependent on its initial spin period Pq, 
which is generally unknown, as £^tot = 27r^/(-^ — 
where / is the moment of inertia of the neutron star. Pg 
can be analytically calculated for an ideal case, assum- 
ing that k and n are constant in the braking equation 
tl = — fcfi" ([Gaensler fc Slanc" 2006). Using the standard 
parameters of the pulsar PSR B1509— 58 (period of 150 
ms, period der ivative of 1.5 x 10~^^ ss ~^ and braking in- 
dex n = 2.84 ([Livingstone et al]|2005| )). we obtain Pq = 
16 ms and Etot = 7.5 x lO'^^ ergs. Knowing that the fit 
requires a total injected energy (i.e. integrated electron 
energies above 1 GeV) of We = 3.0 x 10''^ ergs, 4% of 
i?tot should be converted into the current kinetic energy 
of electrons. 

The equipartition magnetic field strength can be esti- 
mated as Bcq = ^J%nWe/V , where ^ is a volume of the 
emission region. Assuming a spherical region of radius 
r, we obtain Pcq — 22(10pc/r)^/^ /iG. In the case of the 
considered source where r ^ 10 pc, the PW N is particle 
dominated as sugg ested bv [ciTevalieil ([2004D . The mea- 
surement of the extension of the HE 7-ray emission is 
of particular interest to better estimate the boundary of 
the PWN. 

We also consider a 7r° decay model, assuming a pro- 
ton spectrum described by a power-law with a cut-off 
to fit the data points, though there are few theoreti- 
cal indications supporting an injection of such hard pro- 
tons. We obtain an index of protons of 1.9 with a cut-off 
energy of of 60 TeV, which yields the accumulated en- 
ergy of protons above 1 GeV of 1.2 x 10^^(1.0 cvo?)ln 
ergs, where n is a number density of target nuclei. This 
scenario is highly disfavoured from the energetics: even 
wi th a very high density of ^ 10 cm~'^ as mentioned 
by iDubner et al] (|2002D for the northwest limb of MSH 
15—52, the energy required would significantly exceed 
the total energy that the pulsar can supply to its nebula 
[Etot = 7.5 X 10^9 ergs). 

6. CONCLUSION 

We report the detection of pulsed high energy 7-rays 
from PSR B1509-58 below 1 GeV and extended 
emission from its PWN in MSH 15-52 up to 100 GeV 
using 1 year of survey data with Permi-LAT. The 
LAT light curve of PSR B1509-58 above 30 MeV 
presents two peaks. The 7-ray pulse located at phase 
0.33 ± 0.02 is coincident with the main peak observed 
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Fig. 6. — Spectral energy distribution of th e pu lsar wind nebula powered by PSR B1509— 58, from radio to very high energy 7- 
rays. Predicted spectra as discussed in Section 15.21 are overlaid. The total IC spectrum is shown with a solid line while thinner lines 
denote the individual IC components: CMB (dotted), infrared (dot-dashed) and optical (dashed). The dot-dot-dot-dashed line indicates 
the corresponding sy nchrotron emis s ion. A hadronic 7-ray mod el is also overlaid by a thick gray curve. Obser vational data point s: 
ATCA (open circles. ICaensler et al] l|T999h . IGaenslcr ct ah (20021)), BeppoSAX/MECS (magenta thin line region, Mi neo et aTl (j^poT)), 
BeppoSAX/PDS (light green open squares, iFbro t ct al. (2006) and the upper limit is 1 cr), INTE GRAL /IB IS (dark blue open circles, 
IForot et al; (200ff) and the upper limit is 1 ct), RXTE/PCA+HEXTE (dark green dashed line region, 'Marsdcn ct alj UMS)), COMPTEL 
(open diamonds, KuiDcr ct al. (1999)), EGRET (open crosses, [RuiDcr ct al. (1999)), Fermi-LAT (this paper, filled circles), H.E.S.S. (open 
stars, [Aharo nian ct al. (2005)), CANGAROO-III (open triangles, Nakamori ct al. (2008)). Note that COMPTEL and EGRET points 
present the DC emission from this region, which is dominated by the central pulsar for this energy range and could be contaminated by 
the nearby 7-ray pulsar PSR J1509-5850 recently identified by the LAT tAbdo et al.ii2010ai ). 



in X- and soft 7-rays (jUlmer et al.l 119931 : iKuiper et alJ 
Il999| ). The second peak detected at phase 0.96 ± 0.01 
may correspond to the marginal detection reported by 
IKuiper et al] ()1999f ) in the 10 - 30 MeV energy range 
covered by EGRET. The high-altitude emission models 
have problems explaining the peak separation measured 
with Fermi and the radio peak lagging the first 7-peak. 
A confirmatio n of a possible precurs or in the 1351 MHz 
radio profile (jCrawford et al.l 120011) might change this 
interpretation. The 2 a upper limits derived below 
1 GeV confirm the spectral break of PSR B1509-58 in 
the 10-30 MeV energy range. More data are needed 
to measure the pulsar spectrum in the LAT energy 
range. Such estimates can help constrain the shape 
of the spectrum from X- to 7-rays and disentangle 
between the emission models in the magnetos phere of 
PSR B1509-58. Both the high-altitude in odels ijRomanil 
119961: IVenter. Harding fc Guillemotl l2009D and the polar 
cap model ( Harding et al.lll997t ) can accommodate the 
severe spectral break and the low Fermi limits. Details 
such as the break energy and/or X-ray spectral shape 
are not exactly met. The extended 7-ray emission 
observed by the LAT above 1 GeV is spatially coinci- 
dent with the PWN powered by PSR B1509-58. Its 
morphology is well modeled by a uniform disk or a 
Gaussian distribution. The LAT spectrum of the PWN 



above 1 GeV is well described by a power-law with a 
photon index of (1.57 ± 0.17 ± 0.13) and a flux above 
1 GeV of (2.91 ± 0.79 ± 1.35) xlO'^ cm^ s~^. LAT 
analyses of the PWN in MSH 15—52 bring new elements 
to the discussion on the emission models responsible for 
the high to very high emission from this source. The 
hadronic 7-ray scenario is highly disfavoured by the 
new LAT obs ervations, as suggested b y previous TeV 
obser vations ([Aharonian et ahl 120051 : iNakamori et al.l 
I2OO8I ). The multi- wavelength spectrum can be explained 
by synchrotron and inverse Compton processes, as- 
suming a broken power-law spectrum for the electrons. 
The spectral break, constrained by multi- wavelength 
observations, is likely due to an intrinsic break of 
electrons injected from the pulsar wind. About 4 % of 
the pulsar's loss of rotational energy would be required 
to power the 7-rays detected by the LAT, well in the 
range observed for other pulsar wind nebulae. More 
data are required to estimate the detailed morphology 
of the high energy 7-ray emission and better constrain 
the spectral break in the inverse Compton component. 
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